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Abstract 



Vacuum oscillations are considered for the combined solar neutrino obser- 
vations, including the Kamiokande II spectrum data and incorporating theo- 
retical uncertainties and their correlations. Despite the conceptual difficulty of 
the fine tuning between the neutrino parameters and the Sun-Earth distance, 
2-flavor vacuum oscillations provide phenomenologically acceptable solutions. 
There are allowed regions at 99% C.L. for Am^ = (0.45 - 1.2) x 10"^° eV^ 
and sin^ 20 = 0.6 — 1; the best fit solution is x^/d-f. = 19.2/16, which is 
acceptable at 16% C.L. Oscillations for sterile neutrinos are, however, ex- 
cluded by the averaged data at 99.4% C.L. The vacuum oscillation solutions 
predict characteristic energy spectrum distortions and seasonal variations in 
Sudbury Neutrino Observatory, Super-Kamiokande, and BOREXINO. Those 
predictions are given in detail, emphasizing that the vacuum solutions are 
distinguishable from the MSW solutions. 
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Typeset using REVT^ 



I. INTRODUCTION 



The solar neutrino experiments of Homestake Kamiokande and the galhum 

experiments of SAGE and GALLEX show significant deficits of the solar neutrino 
flux when compared to the standard solar model predictions J^jQ; cis summarized in Ta- 
ble |. Among many theoretical propositions to resolve the discrepancy, the astrophysical 
solutions in general are strongly disfavored by the observations P|-P^. Assuming the stan- 
dard properties of neutrinos, the SSMs are excluded 0. Moreover a model independent 
analysis concluded that the lower observed Homestake rate relative to the Kamiokande rate 
is essentially incompatible with all astrophysical solutions ||l3l . On the other hand, the neu- 



trino oscillation hypothesis — the simplest particle physics solution requiring the minimal 
extension of the standard model of particle physics — fits all observations; the Mikheyev- 
Smirnov-Wolfenstein (MSW) fl^ solutions with Am^ ~ 10~^ eV ^ gives a perfect description 
of the data, and the obtained mass range is consistent with the general expectation of the 
grand unified theories with the seesaw mechanism (for recent MSW analyses, see Ref. 



and references therein). Vacuum oscillations [|T6[ also offer phenomenologically acceptable 



solutions. These, however, require a fine tuning between the neutrino parameters and the 
Sun-Earth distance and suffer a conceptual drawback. 

In this paper the combined solar neutrino observations including the Kamiokande II 
spectrum data are investigated for vacuum oscillations (for other recent analysis, see 
Ref. |ll7|-p^). Q Up-to-date experimental data are used. The theoretical uncertainties and 



their correlations are fully incorporated. It is stressed that, since some of the flux uncertain- 



^ Vacuum oscillations were examined recently by Krastev and Petcov for the averaged data of 
the experiments including the theoretical uncertainties, and for the time-divided data without the 



theoretical uncertainties |21|. Here vacuum oscillations are studied without the time-divided data, 
but including the Kamiokande II spectrum data and incorporating the theoretical uncertainties, 
and the consequences for the future experiments are discussed. 
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ties are comparable to the experimental uncertainties, and also they are strongly correlated 
experiment to experiment and flux to flux, a careful treatment of the theoretical uncertain- 
ties is necessary. The omission of the theoretical uncertainties and their correlations can 
lead to a nontrivial change of the allowed parameter space |]T5|. From the parameter space 



obtained from the existing experiments, one can give robust predictions for the next gener- 
ation experiments of Sudbury Neutrino Observatory, Super-Kamiokande, BOREXINO, and 
ICARUS, which will start operating in the mid and late 1990's. In particular the global 
solutions predict characteristic spectrum distortions in SNO and Super-Kamiokande and 
drastic seasonal variations in BOREXINO. Those predictions are explicitly given, and it is 
emphasized that those measurements distinguish the vacuum oscillation solutions from the 
MSW solutions. 

II. GLOBAL ANALYSIS FOR VACUUM OSCILLATIONS 

In 2-state oscillations in vacuum, the survival probability of the electron neutrinos pro- 
duced in the Sun and measured in a detector is given by 

P^^^^^{E,r,R) = 1 -sin2 2^sin2[(/?-r)AmV4^], (1) 

where Am^ and 6 are the squared mass difference and the mixing angle between electron 
neutrinos and the other species of neutrinos to which i/g's oscillate; E is the neutrino energy 
and r and R are the neutrino production location and the detector location measured from 
the center of the Sun. As the reference SSM, the latest model including the helium diffusion 
effect by Bahcall and Pinsonneault is used 0. In constraining the parameter space from 
the observations, the survival probability Pu^-,ue{E,r, R) is convoluted with the neutrino 
production proflle functions [0 and integrated over r, and is also averaged over the Sun- 
Earth distance R when compared to the time-averaged rate of the experiments. For the 
continuous spectrum flux components the survival probability is also integrated over E. 
For the radiochemical detectors the detector cross sections are taken from Ref. 123112^ and 
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included; for Kamiokande the detector cross section, the detector resolution, and the detector 
efficiency are included 0. When the Kamiokande II spectrum data are incorporated, the 
electron spectrum are divided into 14 bins according to the published data 0. 

Since the theoretical uncertainty of the ®B flux (14%) in the Bahcall-Pinsonneault SSM is 
comparable to the experimental uncertainties of Homestake (11%) and Kamiokande (14%), 
it is important to incorporate the theoretical uncertainties of the neutrino fluxes in the 
analysis. The flux uncertainties are strongly correlated experiment to experiment, and also 
flux to flux. The SSM uncertainties are parametrized with the uncertainties of the central 
temperature and the uncertainties of the relevant nuclear reaction cross sections, following 



the prescription described in Ref. |T^. The uncertainties from the detector cross sections 
are also included and are correlated between the two gallium experiments. 

First the analysis has been carried out for Homestake, SAGE, GALLEX, and the averaged 
rate of Kamiokande II and III without the spectrum data. The two gallium experiments are 
treated as independent data points. When we confront observations to a theory, there are 
two statistical questions to ask: 

1. Given the observations, is the theory acceptable or not? 

2. Once we decide that the theory is acceptable, what is the allowed parameter space in 
the theory at a certain confidence level? 

In the present case, the ffist question is quantified by evaluating the goodness-of-fit (GOF) 
for the minimum in the Am^ — sin^ 26 space. The global minimum when the data are 
fit to the oscillation parameters is 5.0 for 2 degrees of freedom ( 2 d.f. = 4 experiments — 2 
parameters). The probability of obtaining the minimum equal to or larger than 5.0 for 2 
d.f. by chance is 8%. That is, the best fit solution for vacuum oscillations is allowed at the 
8% confidence level (C.L.), which is not an excellent fit, but acceptable. There are seven 
local minima in the region /S.m? = (0.35 — 1.3) x 10^^° eV^ and sin^2^ = 0.6 — 1.0; their 
oscillation parameters, x^? and the GOF values are summarized in Table 
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The data have also been fit to Am'^, sin^ 26', and the mean Sun-Earth distance simul- 
taneously. For the best fit solution (Solution D in Table y), the mean Sun- Earth distance 
has to be fine-tuned within 5% to describe the existing observations at 90% CL. Unless we 
invoke some kind of anthropic principle, there is no apparent logical connections between 
the oscillation parameters and the Sun-Earth distance, and this fine tuning is considered as 
a conceptual drawback in vacuum oscillations. 

Once vacuum oscillations are accepted as a phenomenologically viable theory for the 
solar neutrino observations, the CL. regions are obtained by evaluating at each point in 
the Am^ — sin^ 29 space. The allowed regions are defined by 

X'(sin2 Am^) < Xrain + Ax', (2) 

where Xmin is the global minimum, and Ax^ = 4.6, 6.0, and 9.2 for 90, 95, and 99% 
C.L., respectively. The constraints obtained from each Homestake, Kamiokande, and the 
combined gallium experiment at 95% CL. are shown in Fig, |1|. 

When the parameters are fit with the combined observations, the allowed regions are 
essentially controlled by the Homestake data and limited to Am^ = (0.35 — 1.3) x 10"^'^ eV ^ 
and sin^ 2^ = 0.6 — 1.0. The regions allowed by the combined observations at 90, 95, and 
99% CL. are shown in Fig. S; the allowed regions are essentially in agreement with the 



recent analysis by Krastev and Petcov |2l|. The slight differences are mainly due to the 
difference in the input data. 

Since the z/g survival probability in this parameter space (i.e., the Sun-Earth distance ~ 
the oscillation length) are often energy dependent, the analysis has been also carried out 
with the additional Kamiokande II energy spectrum data. The 14 data points obtained 



^This definition of the allowed parameter space assumes a gaussian distribution of the probability 
density around the global minimum. That is, the allowed regions should have elliptical shape, 
and this is only a crude approximation in the present case. Therefore the allowed regions shown 
in Fig. H, |2|, and |^ should be taken as qualitative displays of confidence levels. 
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in the Kamiokande II experiment are used 0; the theoretical uncertainties as well as the 
experimental systematic uncertainties (0.06 x SSM) are correlated among the 14 bins. No 
spectrum data are available from Kamiokande III and the average value is used as an inde- 
pendent data point. The SSM flux uncertainties are also correlated between Kamiokande 
II and III. Q The expected spectrum distortions for three local minima (B, D, and F) for 
Kamiokande along with the actual Kamiokande II data are shown in Fig. ^ The global 
minimum of the combined analysis is 19.2 for 16 d.f. ( = 18 data points - 2 parameters), 
which is acceptable at 84% CL. The neutrino parameters and the GOF for all local 



minima in this region are listed in Table III. The allowed parameter regions are displayed 



in Fig. |. The parameter space allowed at 99% CL. is for 

Am^ = (0.45 - 1.2) X 10"^° eV ^ and sin^ Id = 0.6 - 1.0, (3) 

and remain essentially unchanged by including the Kamiokande II spectrum information. 

Vacuum oscillations can also be considered between electron neutrinos and sterile neu- 
trinos |E5|,p!7lp!9|,pTl . In this case there are no neutral current events from the transformed 



neutrinos in Kamiokande, which would be present and contribute about 15% of the total 
signal in the case of oscillations into or z/^. This absence of the neutral current contribu- 
tions in Kamiokande makes the difference between the Kamiokande rate and the Homestake 
rate effectively wider, and the fits become considerably worse. In the joint analysis with- 
out the Kamiokande spectrum data, the global minimum is 10.1 for 2 d.f., which is an 
extremely poor fit. The GOF is 0.63%, and the solution is excluded at 99.4% C.L. That is, 
including all theoretical and experimental uncertainties, the existing experiments are essen- 
tially incompatible with vacuum oscillation for sterile neutrinos. Q When the Kamiokande II 



^ The possible correlations of the systematic uncertainties in Kamiokande II and Kamiokande III 
were ignored. 

^ The same conclusion was obtained by Krastev and Petcov but with less statistical signifi- 
cance. The difference is entirely due to the difference in the input data. 
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spectrum data are included, the GOF improves (4.8%), but this is an artifact of increasing 
the number of degrees of freedom by adding the Kamiokande II 14 bins, and does not mean 
that the sterile oscillations are favored by the Kamiokande spectrum data; the incompati- 
bility of the data and the sterile oscillations remains essentially unchanged. The oscillation 
parameters and GOF for each local minimum is listed in Table for the averaged rate, 
and in Table for the analysis including the Kamiokande II spectrum data. 

III. PREDICTIONS FOR THE FUTURE SOLAR NEUTRINO EXPERIMENTS 

The oscillation hypothesis for the solar neutrinos, either the MSW mechanism or the 
vacuum oscillations, will be verifiable in the next-generation solar neutrino experiments 
of SNO Super-Kamiokande BOREXINO and ICARUS Moreover those 



experiments should be able to distinguish the vacuum solutions and the MSW solutions, and 
determine the oscillation parameters accurately. The key measurements in those experiments 
are: 

1. The ratio of charged and neutral current events (CC/NC) from the neutrinos in 
SNO. 

2. The energy spectrum measurements in SNO and Super-Kamiokande. 

3. The measurements of time variations of the ^B neutrinos in SNO and Super- 
Kamiokande, and the measurements of time variations of the ''Be neutrinos in BOREX- 
INO. 

The CC/NC measurement in SNO is the gold plated observable in confirming or falsifying 
the oscillation hypothesis. The comparison of the observed value (CC/NC)sno and the SSM 
value (CC/NC)ssM is free from the SSM flux uncertainties, and a depletion of the observed 
charged current signifies the oscillations between and or z/^. Assuming the Bahcall- 
Pinsonneault SSM, the currently allowed parameter space at 95% C.L. predicts 
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(CC/NC)s^o 



0.15-0.5. (4) 



(CC/NC)33^ 

There are two caveats, however. The CC/NC measurement by itself does not distinguish the 
type of oscillations (e.g., whether the MSW effect or vacuum oscillations). Also (CC/NC)sno 
is the same as (CC/NC)ssm if the oscillations are for sterile neutrinos. The measurements 
of the energy spectrum and of time variations are important in determining those detailed 
features of neutrino oscillations. In the following, the predictions for spectrum distortions 
and the time variations are discussed separately, but those two effects are often strongly 
correlated, and the real data in the future experiments should be analyzed accordingly. 

Although the existing Kamiokande energy spectrum data do not significantly constrain 
the parameter space, the next generation experiments with high statistics and with lower 
energy thresholds (~ 5 MeV) are quite sensitive to the spectrum distortions caused by 
oscillations. In Fig. |^, the expected spectrum shapes are displayed for three vacuum solutions 
(B, D, and F defined in Table |T|) along with the spectra expected from the MSW solutions 
obtained by the global analysis; the MSW large-angle solution has little distortion and is 
identical to the astrophysical solutions. The spectra are normalized above the proposed 
detector threshold (5 MeV). The charged current cross section and the detector 

resolutions |^ are also included. The error bars correspond to the statistical uncertainties 
for a total of 6,000 events (equivalent to 2 year operation). The vacuum oscillations are 
clearly distinguishable from the MSW solutions and the astrophysical solutions. The similar 
predictions for Super-Kamiokande are displayed in Fig. ^ with error bars corresponding to 
a total of 16,000 events (2 year operation). 



^In the SNO spectrum analysis, it is important to include the effect of the continuous spectrum 
of the out-going electrons and also the detector resolution: the approximation of the out-going 
electron spectra from the monoenergetic neutrinos with a delta function (i.e., "observed electron 
energy" = "neutrino energy" - 1.44 MeV) significantly overestimates the detector sensitivity for 
spectrum distortions and is inappropriate for the spectrum shape analysis. 
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The survival probability [Eqn. (|T])] is also sensitive to the Sun-Earth distance R and 
the solutions of the combined observation predict seasonal variations of the solar neutrino 
signal due to the eccentricity of the Earth orbit. Especially drastic time variations are 
expected for the ''Be neutrino signals since they are mono-energetic (therefore not smeared 
out by the energy integral), and also the oscillation length is shorter than for most of the 

neutrinos [|l^,|l^,^]. For the existing experiments the expected time variations for the 
global solutions (B, D, and F) are displayed in Fig. 0: the error bars are the statistical 
uncertainties when (overly optimistic) large number of events are assumed: 500, 1,000, and 
1,000 events for Kamiokande, Homestake, and the gallium experiment, while the currently 
accumulated signals are ~ 200, ~ 500, and ~ 50 events, respectively. Obviously there is 
little chance of detecting the seasonal variations by the currently operating experiments. 

The time variations expected in SNO, Super-Kamiokande, and BOREXINO are shown 
in Fig. P: the error bars are for the total number of events 12,000, 32,000, and 12,000, 
respectively (each corresponds to 4 years of operation). The seasonal variations should 
be observable in these experiments. In particular the time variations are most drastic in 
BOREXINO because of its sensitivity for the ''Be neutrinos. Those measurements are un- 
mistakable and should provide strong constraints on the oscillation parameters. 

IV. CONCLUSION 

The Homestake, SAGE, GALLEX and Kamiokande data including the Kamiokande II 
spectrum data are analyzed for neutrino oscillations in vacuum. There are phenomenologi- 
cally acceptable solutions in Aw? = (0.45 — 1.2) x 10"^° eV ^ and sin^ 26 = 0.6 — 1.0, although 
these solutions suffer from the fine tuning between the oscillation parameters and the Sun- 
Earth orbit. Oscillations for sterile neutrinos are strongly disfavored by the observations 
and are excluded at 99.4% C.L. The acceptable solutions for the flavor oscillations predict 
characteristic energy spectrum distortions and seasonal variations in the next generation so- 
lar neutrino experiments. The measurements of these signals should distinguish the vacuum 
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oscillations from the astrophysical solutions as well as from the MSW solutions. 
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TABLES 

TABLE L The standard solar model predictions of Bahcall and Pinsonneault (BP SSM) 
and of Turck-Chieze and Lopes Q (TCL SSM). The Bahcall-Pinsonneault model includes the 
helium diffusion effect, while the Turck-Chieze-Lopes model does not. Also shown are the results 
of the solar neutrino experiments. The gallium experiment is the combined result of SAGE and 
GALLEX I and IL 

BP SSM TCL SSM Experiments 

Kamiokande ^ 1 ± 0.14 0.77 ± 0.19 0.51 ± 0.07 BP SSM 

Homestake ^ (SNU) 8 ± 1 6.4 ± 1.4 2.34 ± 0.26 (0.29 ± 0.03 BP SSM) 

Ga experiment (SNU) 131.5 tl 122.5 ±7 81 ± 13 (0.62 ± 0.10 BP SSM) 

SAGE " (SNU) 70 ± 22 (0.53 ± 0.16 BP SSM) 

GALLEX ^ (SNU) 87 ± 16 (0.66 ± 0.12 BP SSM) 

^The result of the combined data of 1040 days of Kamiokande II [0.47 it 0.05 (stat) it 0.06 (sys) 
BP SSM] § and 514.5 days of Kamiokande III [0.57 ± 0.06 (stat) ± 0.06 (sys) BP SSM] g. 
•^The result of Run 18 to 122 is 2.34 ± 0.16 (stat) ± 0.21 (sys) SNU §. 

=The prehminary result of SAGE I (from January, 1990 through May, 1992) is 70 ± 19 (stat) ± 10 
(sys) SNU [|. 

'^The combined result of GALLEX I and II (including 21 runs through April, 1993) is 87 ± 14 
(stat) ± 7 (sys) SNU @. 
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TABLE II. The solutions for vacuum oscillations when the combined result of Homestake, 
Kamiokande, SAGE, and GALLEX is used. The SAGE and GALLEX data are fit separately. 
The Kamiokande II spectrum data are not included. A - G correspond to seven local ^ minima 
within the allowed regions at 99% C.L. The goodness-of-fit (GOF) is the probability of obtaining 
by chance a minimum value equal to or larger than the value we actually obtained. 



Am^ (eV2) sin2 2Q ^ (2 d.f.) GOF (%) 

A 1.2 X 10-1° 1.00 10.9 0.42 

B 1.0 X 10^1° 0.96 5.7 5.9 

C 9.0 X 10-11 0.84 5.1 7.9 

D 7.9 X 10-11 0.78 5.0 8.0 

E 6.3 X 10-11 0.81 5.2 7.6 

F 5.2 X 10-11 0.99 7.5 2.3 

G 3.8 X 10-11 1.00 14.5 0.071 



TABLE III. Same as Table || except that the Kamiokande II spectrum data are included. 



Am2 (eV^) sin^ 20 ^ (16 d.f.) GOF (%) 

A 1.2 X 10-1° 1.00 25.4 3.1 

B 1.0 X 10-1° 0.91 20.5 12 

C 9.1 X 10-11 0.82 19.2 16 

D 7.8 X 10-11 0.78 19.2 16 

E 6.5 X 10-11 0.80 20.4 12 

F 5.2 X 10-11 0.99 23.7 4.9 

G 3.8 X 10-11 1.00 29.9 0.79 
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TABLE IV. The solutions for vacuum oscillations for sterile neutrinos when the combined 
result of Homestake, Kamiokande, SAGE, and GALLEX is used. The SAGE and GALLEX data 
are fit separately. The Kamiokande II spectrum data are not included. A - G correspond to seven 
local minima within the region. 



Am2 (eV2) sin2 26 (2 d.f.) GOF (%) 

A 1.2 X 10-1° 1.00 13.5 0.12 

B 1.1 X 10-1° 0.88 10.1 0.63 

C 9.1 X 10-11 0.79 11.0 0.42 

D 7.8 X 10-11 0.75 11.0 0.41 

E 6.3 X 10-11 0.78 10.4 0.56 

F 5.0 X 10-11 0.93 12.7 0.18 

G 3.8 X 10-11 1.00 16.5 0.026 



TABLE V. Same as Table IV except that the Kamiokande II spectrum data are included. 



Am^ (eV^ 



sin^ 29 



(16 d.f.) 



GOF (%) 



A 
B 
C 
D 
E 
F 
G 



1.2 X 10-1° 
1.1 X 10-1° 
9.1 X 10-11 
7.8 X 10-11 

6.3 X 10-11 
5.0 X 10-11 



3.8 X 10 



-11 



0.98 
0.84 
0.78 
0.74 
0.77 
0.91 
1.00 



30.1 
27.0 
27.6 
27.1 
26.4 
29.3 
32.3 



1.7 
4.2 
3.5 
4.1 
4.8 
2.2 
0.92 
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FIGURES 

FIG. 1. The allowed parameter space at 95% C.L. for the Homestake, Kamiokande and 
combined gallium experiments of SAGE and GALLEX including the theoretical uncertainties and 
their correlations. The Kamiokande II spectrum data are not included. 

FIG. 2. The parameter space allowed by the combined observation of the Homestake, 
Kamiokande, and the combined gallium experiments of SAGE and GALLEX, including the theo- 
retical uncertainties and their correlations. The Kamiokande II spectrum data are not included. 

FIG. 3. The spectrum shape divide by the SSM spectrum expected in Kamiokande. The 



shape is shown for the combined solutions B, D, and F, which are listed in Table III. Also shown 
are the Kamiokande II spectrum data |^] 

FIG. 4. The parameter space allowed by the combined observation of the Homestake, 
Kamiokande, and the combined gallium experiments of SAGE and GALLEX, including the theo- 
retical uncertainties and their correlations. The Kamiokande II spectrum data are included. 

FIG. 5. The electron spectrum shape expected in SNO charged current reaction for the vac- 



uum solutions B, D, and F (defined in Table III). Also shown are the spectra expected for the 



MSW solutions obtained from the combined observations. The spectrum shape of the astrophys- 
ical solutions (or no oscillations) is identical to the MSW large-angle spectrum. The spectra are 
normalized above the threshold (5 MeV). The error bars correspond to the statistical uncertainties 
for a total of 6,000 events (2 year operation). 
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FIG. 6. The electron spectrum shape expected in Super-Kamiokande for the vacuum solutions 



B, D, and F (defined in Table III). Also shown are the spectra expected for the MSW solutions 
obtained from the combined observations. The spectrum shape of the astrophysical solutions (or 
no oscillations) is identical to the MSW large-angle spectrum. The spectra are normalized above 
the threshold (5 MeV). The error bars correspond to the statistical uncertainties for a total of 
16,000 events (2 year operation). 



FIG. 7. The seasonal variation of the signal expected in the existing experiments for the 
vacuum solutions B, D, and F (see Table III). The error bars corresponds to (overly) optimistic 
statistical errors. The accumulated events in the existing experiments are ~ 200, ~ 500, and ~ 50 
for Kamiokande, Homestake, and the combined gallium experiments. No significant constraint is 
expected from the existing experiments. 



FIG. 8. The seasonal variation of the signal expected in the future experiments of SNO 
(charged-current reaction), Super-Kamiokande, and BOREXINO for the vacuum solutions B, D, 
and F (see Table |l|). The error bars correspond to the total events for 4 year operation. Those 
measurements should be sensitive enough to distinguish the vacuum oscillations from the MSW 
solutions and constrain the parameter space. 
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